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Abstract 

The ab initio self-interaction-corrected local-spin-density approximation is used to study the 
electronic structure of both stoichiometric and non-stoichiometric nickelates. From total energy 
considerations it emerges that, in their ground-state, both LiNi02 and NaNi02 are insulators, with 
the Ni ion in the Ni3+ low spin state (i^g^g) configuration. It is established that a substitution of a 
number of Li/Na atoms by divalent impurities drives an equivalent number of Ni ions in the Ni02 
layers from the JT-active trivalent low-spin state to the JT-inactive divalent state. We describe 
how the observed considerable differences between LiNi02 and NaNi02 can be explained through 
the creation of Ni^+ impurities in LiNi02. The indications are that the random distribution of 
the Ni'^"'" impurities might be responsible for the destruction of the long-range orbital ordering in 
LiNi02. 
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The metal oxide compounds LiM02 (M=Ni, Co, Mn) are of considerable interest due 
to their potential application as a cathode material in lithium-ion rechargeable batteries^. 
In particular, LiNi02 and its isoelectronic counterpart NaNi02 have been the subject of 
numerous model and first-principles calculations focussing on the possible charge, spin and 
orbital degrees of freedom, and their impact on the electronic and magnetic properties of 
these compounds. In most model calculations, a JT-active Ni^+ ion with ^29^9 ground-state 
configuration has been assumed, in agreement with experimental evidence — , although a 
Ni^"'"((i^) ion has also been suggested^. With the Ni'^+ configuration, model calculations 
have been able to reproduce the experimentally observed orbital ordering in NaNi02 and for 
LiNi02, a number of models have been proposed to explain the absence of long range mag- 
netic order (see for example Reynaud et al^ and references therein) . However, no universally 
accepted model has emerged that can simultaneously explain magnetic and orbital ordering 
in NaNi02 and the absence of the very same in LiNi02. First-principles band structure cal- 
culations based on the local spin density (LSD) approximation to density functional theory 
(DFT) have not even been able to yield the expected insulating ground state^. Assuming 
a low-spin ground-state configuration for the Ni^+-ion, electronic structure calculations, us- 
ing the LDA+U approach, result in the opening of an insulating gap in both LiNiOj^ and 
NaNiOa^. Although very different in nature from model approaches, the LDA-I-U method 
also relies on the introduction of the Hubbard U to account for the strong correlations. As 
there is no straightforward approach for determining the U from first-principles, it is usually 
treated parameter in the calculations. 

The motivation for the present letter is to apply the self-interaction corrected (SIC) LSD 
method and seek first-principles, parameter free, understanding of the ground state spin 
and valence configurations of nickelates. Based purely on the total energy considerations, 
we show that Ni^+ ground-state configuration in the stoichiometric compounds. 

Furthermore we predict valency changes upon substitution of Li for Ni with potentially far 
reaching consequences for the magnetic and orbital odering. 

The SIC-LSD method^ is an ah initio approach that corrects for an unphysical self- 
interaction contained in the LSD total energy functional^". While, for extended band states, 
the resultant error in the LSD energy is generally insignificant, it may be considerable for 
atomic-like localized states. In the SIC-LSD method both localized and delocalized states 
are expanded in the same set of basis functions, and are thus treated on an equal footing. 
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Different localized/delocalized configurations are realized by assuming different numbers of 
localized states - here (i-states on Ni-atom sites. Since the different localization scenarios 
constitute distinct local minima of the same energy functional, E^^'-^ , their total energies 
may be compared and the global energy minimum then defines the ground state total energy 
and the valence configuration of the Ni-ion. This latter is defined as the integer number of 
electrons available for band formation, namely 



val 



Z — Nrnrp. 



SIC, 



(1) 



where Z is the atomic number (28 for Ni), Ncore is the number of core (and semicore) electrons 
(18 for Ni), and Nsic is the number of localized, i.e. self-interaction corrected, states. 

LiNi02 crystallizes in a trigonal structure {R3m (166) space group) that consists of layers 
of Ni02 slabs built from edge sharing NiOe octahedra and separated from each other by a 
layer of Li cation o^^i^^ (see Fig. 1). This crystal structure, which can be derived from 



.a. 




FIG. 1: LiNi02 crystal structure, consisting of alternating layers of O (blue), Li (red), and Ni 
(green) atoms. 



the rock salt structure of NiO by replacing every second Ni (111) plane by a Li plane, is 
favorable to fast-ion intercalation/removal of Li ions during a battery discharge/charge cycle. 



NaNi02 is isoelectronic to LiNi02 and, at high temperatures, also crystalhzes in the R3m 
structure. However, below T=480 K, NaNi02 undergoes a crystal structure distortion from 
rhombohedral to monoclinic (C2/m (12) spacegroup) that is driven by a cooperative Jahn- 
Teller ordering in the Ni-ion layer—. In LiNi02, a local distortion of the NiOg tetrahedra 
has been observeciii and associated with the expected Jahn- Teller (JT) activity of the Ni 
ions in an octahedral environment. However, unlike NaNi02, the crystal structure shows 
no indication of collective JT-distortion. It has been suggested, that this is due to the 
formation of three sublattices with non-coUinear JT-axes, and the strain fields generated by 
the triangular symmetry of the Ni02 slabs prevent the local orbitral ordering from developing 
into long-range ordering—. The difference in magnetic behaviour is equally puzzling in that 
the observed magnetic order of NaNi02, ferromagnetic intra-(Ni)layer coupling with anti- 
ferromagnetic inter-(Ni)layer couplin g^^i^^ , does not seem to occur in LiNi02. 

Assuming the undistorted R3m crystal structure for LiNi02, we have calculated the total 
energies for five Ni valences - Ni^"*" through Ni^^ - as well as the completely delocalized 
LSD configuration. To determine the lowest total energy for each Ni valency (i.e. number 
of locahzed d-states), we have considered several possible ways the localized states can be 
permuted among the available d-orbitals. The results are summarized in Fig. |21 Clearly, 
the LSD is the least favorable scenario. Overall, the different valency configurations become 
energetically more favorable as the number of localized rf-electrons increases. The global 
energy minimum is obtained for the trivalent Ni'^^ scenario with seven localized d-electrons 
in the t2ge^ low spin configuration. This is marked by the green line in Fig. El Compared to 
Ni^"*", the Ni^"*" configuration is again less favorable. Based on total energy considerations, 
the SIC-LSD thus validates the choice of charge and spin configuration used in the model and 
LDA+U calculations of LiNi02. The same groundstate configuration is found in a similar 
study of the electronic structure of NaNi02. 

In actuality an indication of the SIC-LSD preference for the low spin ground state of the 
Ni-ion is already inherent in the LSD calculation. As can be seen from the density of states 
(DOS) shown in Fig. IHK, the splitting of ~ 2.0 eV between the majority Ni t2g and eg states 
(crystal field splitting) is some four times larger than the exchange splitting (~ 0.5 eV) 
between majority and minority Ni (i-states. A detailed inspection of this same DOS reveals 
that the states in both spin channels are fully occupied while the majority eg states are 
only partially filled, as the Fermi energy falls in the range of hybridized O 2p and Ni eg 
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FIG. 2: Total energy as a function of Ni d-configuration in LiNi02. A number of scenarios 
are shown for each given configuration, dlydlzdlzdl^2_r2dl2_y2 - diydyzdizdl^^2_j.2d[.2_y2, where 
a locahzed spin up/spin down state is indicated by the j/J, symbol, and a delocalized state is 
indicated by the symbol. 

states. Since the self-interaction is only sizeable for the occupied and fairly localized states, 
one can easily see that only these occupied Ni d-states will benefit from the self-interaction 
correction, thus leading to the global energy minimum of the SIC-LSD energy functional for 
this low-spin state. 

When comparing the SIC-LSD DOS for the ground state configuration (Fig. Eb) to the 
LSD DOS (Fig. EK), one notices that the t2g states at the top of the valence band and parts of 
the spin-up eg band states at the bottom of the conduction band have vanished. The Fermi 
level is now situated above the completely filled valence band and LiNi02 is an insulator, 
in agreement with experiment. Because SIC-LSD is still one-electron ground state theory 
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FIG. 3: Spin-decomposed density of states as obtained for a) LiNi02 with the LSD, b) LiNi02 
with the Ni2+ groundstate configuration, c) LisNi^+Ni^+Os, and d) LisNi^+Ni^+Ni^+Og. The 
green and black hnes represent respectively the Ni d-projected and total DOS. The red in c) and 
d) represents the d-projected DOS of the Ni substituted in the Li layer. The perpendicular dashed 
blue line indicates the position of the Fermi level. 

it does not give accurate removal energies of localized states, due to the unaccounted for 
screening and relaxation effects^^. The localized states are therefore situated unphysically 
low at around -0.8 Ry (not shown). 

Given the similarity of electronic structures of stoichiometric LiNi02 and NaNi02, 
what then is the origin of their very different properties? Although super-stoichiometric 
Lii_a.Nii+j.02, < X < 0.2, was first synthesized by Dyer et alM- in 1953, attempts to 
synthesize the stoichiometric LiNi02 {x = 0) compound have so far proved unsuccess- 
ful. There is now a growing body of evidence that chemical disorder inherent to super- 
stoichiometric LiNi02, but absent in NaNi02, is an important factor contributing to prop- 
erty differences. Indeed, in a recent theoretical paper, Mostovoy and Khomskii^^ raise the 
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possibility that the complex electronic properties of LiNi02 are related to the presence of 
Ni anti-site ions in the Li layers. A number of experimental studies have also come to 
the same conclusio n^^i^" . In super-stoichiometric Li1_j.Ni1_1_2.O2 with excess Ni- atoms pre- 
sumed to occupy the Li-layer, charge balance arguments suggest that misplaced Ni-atoms 
will have a 2+ charge state ^2iSi which will also induce an equivalent number of Ni^^-ions 
in the Ni02 layer. Namely, a valence state that can be captured by the crystallographic 
formula (Lilt^^il^) Li-iayeri^^l'^^iitrc) Ni-iayer02. Despite its apparent complexity, the en- 
ergetic basis of such a ground state can be investigated straightforwardly using the SIC-LSD 
method. 

Starting from a supercell consisting of four LiNi02 formula units (crystal structure R3m) , 
we replace a single Li atom with Ni. In Table I, we compare the total energy for the trivalent 
scenario, i.e. (Lil'^'Ni^^) n^iayeri^il'^) Ni~iayer02 in row 1, to a number of scenarios where an 
increasing number of Ni atoms is treated as divalent Ni^"''((i^). From Fig. we notice that in 
the energetically unfavorable, all Ni^"*" configuration, the substitution of a monovalent Li by a 
trivalent Ni results in additional e^ bands that are partially filled. The gain in hybridization 
energy associated with this d-band formation is however not large enough to compete with 
the possible gain in self-interaction energy that results from localizing the corresponding 
d-states. Consequently, the energetically most favorable scenario is obtained with two Ni 
atoms in the Ni^"^ configuration and three Ni atoms in the Ni'^+ configuration (row 4 in 
Table I), i.e. {Lil'^'Ni^'^) Li_iayer(^^l^^^'^^)Ni-iayer02- The corresponding DOS is shown in 
Fig. inil, where it can be seen that, due to the localization of d-states, a large band gap is 
observed indicating that the compound remains insulating even when off-stoichiometric. 

In the preceding discussion it is important to keep in mind that the valencies obtained 
in SIC-LSD are not ionic in nature. Rather, valency change is the result of a rf-electron 
localization/delocalization transition and is not necessarily indicative of large charge trans- 
fer. Whether a given d-state will localize or take part in d-band formation depends on the 
outcome of the competition between hybridization and self-interaction energies. 

The above interesting behavior of misplaced Ni-atoms notwithstanding, the question still 
remains as to the extent that the off-stoichiometry induced valence change can explain the 
observed differences in electronic and magnetic properties between LiNi02 and NaNi02. 
Unlike Ni^"*" ions, divalent Ni^"*" ions are not JT-active. Consequently, a random distribution 
of these ions in the Ni02 layer could possibly affect the orbital ordering. However, this 
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TABLE I: Total energies for different Ni^+/Ni^+ configurations in LisNisOg, and Na3MgNi408. 
The energies are given relative to the corresponding all trivalent Ni^"^ configuration. The Ni/,j 
notation refers to a Ni substituted for a Li atom. 



LigNisOg 


Ni valency 


Total energy (Ry) 


Ni3+ Ni3+ Ni3+ Ni3+ 

Li 


0.00 


Ni3+ Ni3+ Ni3+ Ni3+ Ni2+ 

Li 


-0.0393 


Ni2+ Ni3+ Ni3+ Ni3+ Ni3+ 


-0.0691 


Ni^rt Ni3+ Ni3+ Ni3+ Ni2+ 


-0.1033 


Ni2+ Ni3+ Ni3+ Ni2+ Ni2+ 


0.0262 


Ni2+ Ni2+ Ni2+ Ni2+ Ni2+ 

Li 


0.1494 


Na3MgNi408 


Ni valency 


Total energy (Ry) 


Ni3+ Ni3+ Ni3+ Ni3+ 


0.00 


Ni3+ Ni3+ Ni3+ Ni2+ 


-0.0136 


Ni3+ p^i3+ p^j2+ j^-2+ 


+0.0658 



issue is clouded by the fact that introduction of Ni ions into the Li layer gives rise to 
additional magnetic interactions between neighboring Ni layers which then need to be taken 
into account on top of the pre-existing intra-layer magnetic interactions. For example, it 
has been suggested that the substitution of some of the Na atoms by Ni atoms in NaNi02 
leads to frustration of the antiferromagnetic stacking order and the destruction of long range 
magnetic order—. Given this complexity, it would clearly be advantageous to separate the 
effect of valency change on the orbital ordering and on the magnetic interactions. In what 
follows we suggest how this may be achieved. 

Specifically, we propose to start from stoichiometric NaNi02, which is both magnetically 
and orbitally ordered, and substitute Mg or Ca ions for some of the Na ions. By replacing 
monovalent Na with divalent Mg/Ca, we expect no additional interlayer magnetic interac- 
tions to occur, but the extra valence electron will, dependent on the relative size of the 
corresponding self-interaction and hybridization energies, either localize on a Ni-site or oc- 
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cupy an conduction band level. Starting from a Na4Ni408 unit cell and replacing a single 
Na atom by Mg, we have studied a number of Ni-ion configurations. As shown in Table I, 
we find that the configuration with three Ni ions in the 3+ charge state, and one Ni ion 
in the 2+ charge state is energetically most favorable. Thus, the substitution of Mg/Ca 
atoms in the Na layer, leads to localization on an equivalent number of Ni-ions, and again 
results in a mixture of JT active and non-active ions in the Ni02 layer. Clearly, investigating 
experimentally how this would affect the local and medium range orbital ordering is then of 
considerable interest. 

A study — of mixed layered oxide phases Nai_2,.Li2.Ni02, where Li is gradually replaced 
by Na, has shown that both orbital-ordering and antiferromagnetism are conserved for < 
X < 0.2, and that the collective Jahn- Teller transition is only suppressed for x > 0.2. This 
result is in line with our picture, as the substitution of Na by isoelectronic Li does not result 
in a change in valency of the Ni ions in the Ni02 layer. If the substitution of Mg or Ca for 
Na brought about the suppression of long range orbital ordering for comparatively smaller 
concentrations, this would be a strong indication that it is indeed the random distribution of 
non JT-active Ni ions in the Ni02 layer that is responsible, rather than the introduction of 
magnetic Ni atoms in the Li layers. On the other hand if the orbital ordering persists up to 
higher Mg/Ca concentrations, it is more likely that the small ionic size of Li"*" is responsible 
for the suppression of macroscopic Jahn- Teller distortion. 

In summary, we have applied the SIC-LSD method to study the electronic structure of 
both stoichiometric and off-stoichiometric nickelates. Based on parameter free total energy 
considerations, we find the ground-states of both LiNi02 and NaNi02 to be insulating, 
with Ni^"*" in the tg^e^ low-spin configuration. Furthermore, the calculations reveal that the 
crystal field effects rather than the exchange interactions determine the low spin ground 
state configuration. With respect to Lii_a;Nii+a;02, we find that for every Ni substituting 
for Li, a Ni-ion turns divalent in the Ni02 layer. The presence of a divalent Ni-ion will 
alter the exchange interaction, and would strongly affect the magnetic properties of these 
compounds. We propose that NaNi02 having had Mg or Ca substituted for some of its Na 
would be ideally suited to experimental study of the effect of Ni^"*" ions on the long range 
orbital ordering. 
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